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Abstract 

We consider the effect of a random longitudinal field on the Ising model in 

a transverse magnetic field. For spatial dimension d > 2, there is at low 

strength of randomness and transverse field, a phase with true long range 

order which is destroyed at higher values of the randomness or transverse 

field. The properties of the quantum phase transition at zero temperature are 

controlled by a fixed point with no quantum fluctuations. This fixed point 

also controls the classical finite temperature phase transition in this model. 

Many critical properties of the quantum transition are therefore identical to 

those of the classical transition. In particular, we argue that the dynamical 

scaling is activated, i.e, the logarithm of the diverging time scale rises as a 

power of the diverging length scale. 
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A number of recent theoretical [|T]-§[ and experimental |0] works have studied the effects 
of randomness on simple quantum statistical models. These studies have been motivated 
by the need to understand in a simple context the interplay of effects related to strong 
randomness, interactions, and quantum fluctuations. In this paper, we study the effect 
of a random longitudinal field applied to the transverse field Ising model. The effects of 
random longitudinal fields on classical Ising models have been studied extensively [§] , and 
are partially well understood. In contrast, very little is known about the corresponding 
quantum problem in realistic dimensions. Most previous studies of this system have been 
limited to mean field theory (expected to be valid above 6 spatial dimensions; see below) 
or to an expansion in e = 6 — d. However, as is well known from the classical problem, 
these results are not expected to directly be of much relevance to physical systems in finite 
dimensions much smaller than 6. Here we will provide a general scaling theory applicable 
in any dimension. 

The model in consideration is defined by the Hamiltonian 

^ = E ^1 - ^ E - E (1) 

i,j i i 

For simplicity, we have assumed that J and h are non-random, though our results should 
apply also for weakly random J and h. Hi is assumed to be random with zero mean and 
variance A^. The ai are Pauli spin matrices. We note that, as usual, this quantum model 
in d dimensions at T = is equivalent to a classical Ising model in d + 1 dimensions in a 
random field with the randomness correlated along one direction. 

It will often be more convenient to consider a coarse-grained continuum field theoretic 
version of the Hamiltonian Eqn. |1|. The continuum action is readily written down as: 
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S = I dr I d'^x-^ 



(2) 



where H{x) is the (coarse-grained) random field and is taken to be Gaussian distributed 
with mean and variance A^. For future convenience, we have introduced a factor t as an 
overall scale for the action. 



Later in this paper, we use the standard Imry-Ma argument to show that the ordered 
phase is unstable to any weak random field for spatial dimension d < 2. In Fig |1], we 
show the expected phase diagram for d > 2 in the temperature (T), h, A space. In what 
follows, we will concentrate primarily on the T = transition from the ordered phase to the 
paramagnetic phase. Our central claim is that the quantum transition is controlled by a fixed 
point with no quantum fiuctuations in any dimension (> 2). However quantum fiuctuations 
are "dangerously irrelevant" at this fixed point, and need to be included to get the correct 
critical behaviour of many physical observables. Such fiuctuationless fixed points have arisen 
in several other recent studies of quantum transitions in random systems ||T|,p|,^|6|,p!0| . For the 
model in consideration here, this claim is completely analogous to the corresponding claim 
on the role of thermal fiuctuations at the finite temperature phase transition in classical 
random field Ising magnets . Indeed, we claim that the same fixed point controls both the 
classical (finite T) and quantum (T = 0) transitions in this model. This enables us to relate 
many of the critical properties of the quantum transition to those of the well-studied (but 
only partially understood) classical transition. In particular, we suggest that the dynamical 
scaling is activated, i.e, the logarithm of the diverging time scale rises as a power of the 
diverging length scale. 



I. WEAK RANDOMNESS 

We first consider the effects of weak random fields on the properties of the pure system. 
These are quite innocuous deep in the paramagnetic phase when the transverse field wins 
over both the randomness and the exchange interaction. So we will turn immediately to the 
ordered phase and the critical point. 



It is easy to see using the Imry-Ma [ITTHT^ argument that for dimension d > 2, the 



ordered phase of the quantum model is stable to a weakly random Hi. Imagine that in the 
absence of randomness, the system is deep in the ordered phase. The stability to a weak 
random field is determined by balancing the energy cost to form large-sized domains with 
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the typical energy gain because of the random Hi. Deep in the ordered phase, the energy 
cost of a domain ~ L'^~^ while the typical energy gain due to the random field ~ L'^/^ 
exactly as in the classical problem. Thus for > 2, for weak enough random fields, it is not 
favourable to form large-sized domains and the system stays ordered. We expect of course 
that strong randomness will destroy the ordered phase in any dimension. Precisely in d = 2, 
as in the classical problem, the system is marginally unstable and there is no ordered phase. 

Now, assume that we are at the critical point of the pure system ind > 2. When are weak 
random fields a relevant perturbation at this critical fixed point? This can be answered as 
follows. Consider the continuum action Eqn. ^j. Averaging over the disorder using replicas 
gives the term ^ J d'^x J dTidT2j2a,b 'Paix, Ti)(f)h{x, T2) . Now under the renormaliation group 
transformation appropriate to the critical fixed point of the pure system x — x' = |, 
r — > r' = ^, and — 0' = 0s 2 . This gives A^' = A^s^^^"*^. Clearly then A is always 
relevant. 



II. THE CRITICAL FIXED POINT 



In this section, we will present the reasoning behind our assumption that the fixed point 
controlling the transition has no quantum fluctuations . 



For classical random field magnets, it is believed [Q that the properties of the finite 
temperature phase transition (which exists for d > 2) are controlled by a zero temperature 



fixed point. The first evidence for this belief came from perturbative studies |[TT| , p!5[| of 
a continuum field theoretic description of the magnet near the critical point. Order by 
order in perturbation theory, it can be shown that the effects of the fluctuations introduced 
by the randomness dominate over the effects of the thermal fluctuations near the critical 
point. This is interpreted in renormalization group language to be a manifestation of a 
zero-temperature fixed point with (dangerously) irrelevant thermal fluctuations. Though 
many of the predictions of this perturbative analysis are believed not to be correct for 
realistic dimensions {d = 2,3), there is a general consensus that the critical fixed point is 
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fluctuationless in all dimensions where there is a transition. This is further supported by 



heuristic renormalization group arguments valid in ci = 2 + e |T^, numerical real space 
renormalization calculations [|T^, and general agreement of some of the phenomenology 
suggested by such fluctuation-less fixed points with experiments in = 3 |T^. 

For the quantum random field magnets that we are considering here, a similar perturba- 
tive analysis, valid in high dimension, of a continuum field theory expected to describe the 
correct T = critical behaviour was undertaken several years ago by two groups |]19|,^ (See 



also Section |rV| ). Again it was found that order by order in perturbation theory, the effects 
of fluctuations introduced by the randomness dominated the effects of quantum fluctuations. 
We take this to be strong evidence that the fixed point is fluctuationless for every dimension 
{d > 2) for the quantum problem as well. 

Now consider the finite temperature phase transition in the quantum model which occurs 
when the ground state is ordered. This transition is, of course, in the same universality class 
as that in the classical models discussed earlier. It is clear that this is controlled by the same 
fixed point that controls the quantum transition. Thus, we have the unusual situation that 
the same fixed point controls both the T = and T ^ transitions. The phase diagram 
and RG flows (for d > 2) are shown in Figure |I|. 



III. GENERAL SCALING HYPOTHESIS FOR 2 < L> < 6 

A number of results follow from the claim that the same fixed point controls both the 
classical and quantum transitions. First, on approaching the transition from the ordered 
phase, the magnetization vanishes with an exponent jSq which is the same as for the classical 
transition. Next consider correlation functions. As in the classical problem, there are two 
different correlation functions that can be defined: The correlation function 

C{x,t;x',t') = < o-2(x, r) >< cr2(x',r') > - < ^^(x, r) > < a;,{x',T') > (3) 

where the angular brackets denote averaging over quantum fluctuations and the overline 
denotes averaging over the randomness. As the randomness is independant of time, the 
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system is translationally invariant in time in every sample. Therefore C{x, r; x', t') is inde- 
pendant of r, r' and we will refer to it as C{x,x') from now on. After averaging over the 
disorder, spatial translational invariance is also restored in the correlation function, and so 
C{x, x') = C{x — x'). This satisfies (near the transition) 

C(x) ~ -^C (l) (4) 

where the correlation length ^ ~ \h — hc\~'^^ ■ The scaling function C and the exponents 
uq, t]q are properties of the fixed point theory and its relevant perturbations (the deviation 
from the critical randomness strength). As these are the same for both the quantum and 
classical transitions, we get the result that the exponents uq, r/g, and the function C are 
identical to their classical counterparts. 
The connected correlation function 

Gc{x, r; x', r ) = < T,(o-^(x, r)a^{x', r')) > - < cr,{x, r) >< a^{x', r') > (5) 

where is the time-ordering symbol. Clearly Gc{x, r; x', r') = Gc{x — x' ,t — t'). 

Note that this correlation function vanishes at the fixed point (as there are no quantum 
fluctuations at the fixed point). Thus to obtain its critical behaviour, we need to keep 
the irrelevant quantum fluctuations. Similarly for the corresponding classical problem, it is 
necessary to keep the irrelevant thermal fluctuations. As these two irrelevant perturbations 
may have quite different effects, in general, we do not expect Gc to scale identically to the 
classical space and time dependant correlation function. 

However it is possible to argue that the static correlation function of the quantum problem 
scales identically to the equal-time correlation function of the classical problem. This was first 
noticed within the e expansion by Boyanovsky and Cardy , but as we show below is more 
generally valid (though other results of the e expansion are not). Consider the susceptibility 
||2T|of the model to an external spatially varying static magnetic field Hext{x)- Clearly for 
both the quantum and classical problems, the scaling of this susceptibility near the transition 
is just determined by minimizing the classical fixed-point Hamiltonian in the presence of the 
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external field. Hence the static (non-local) susceptibilility scales identically for the classical 
and quantum problems. Now at any finite T 7^ in the classical problem, the fluctuation- 
dissipation theorem implies that the equal-time correlation function is proportional to the 
susceptibility. Similarly, the static correlation function of the quantum problem is also 
proportional to its susceptibility. Thus these two correlation functions scale identically. We 
then have the result 

/^-Gc(x,r)~^^(|) (6) 

with 1] and the scaling function Q being identical to those for the connected equal-time 
correlation functions of the classical problem. The equal-time correlation function in the 
quantum problem will however scale differently. 

We now turn to dynamical correlations. For the classical transition. Villain and Fisher 



24 have presented arguments to show that the dynamical scaling is unconventional, with 
the logarithm of the characteristic relaxation times scaling as a power of the length scales. 
Similar unconventional dynamic scaling also occurs in the vicinity of the random quantum 
transitions studied recently all of which are controlled by fluctuation- less fixed points. 

We suggest below that the dynamic scaling is activated at this quantum transition as well. 

The argument for the dynamic scaling closely follows that for the classical case. Consider 
a block of the system of size ~ ^ near the critical system. Ignoring quantum fluctuations, 
the energy landscape as a function of the total magnetization of the block has been argued 
to scale as [^]. Most blocks would thus have a single deep global minimum at some 
non-zero value of the magnetization. The energy of this minimum will differ from those of 
other local mimima by amounts ~ The barriers separating these different local mimima 
also scale as Effects of quantum fluctuations on such blocks should be rather small, and 
do not contribute significantly to the dynamics. However, there would be some rare blocks in 
which there are two minima with an energy separation which is nearly zero. The dynamics 
at long time scales will be dominated by quantum tunnelling between such minima in these 
rare blocks. As the barrier between these minima rises as a power of ^, it is natural that 



the tunneling time tg ~ exp(c^^) with a new exponent. (This quantum tunnehng will 
be important so long as the energy difference between the two minima is roughly less than 
h/tq). Thus the dynamics is activated like in the other random quantum transitions studied 
in Refs 

This therefore motivates the following scaling form for the imaginary part of the g, uo de- 
pendant susceptibility (which is the spectral density for the connected correlator introduced 
above): 

x"(?,^) = r/(ge,!^) (7) 

where to is a microscopic time scale. The exponent k will be related to other exponents 
below. The real part of the susceptibility may be obtained by the Kramers-Kronig relation: 



In In 

Following Pytte and Imry |23|], we write y = ° and z = ■ Then 



TT U' — U 

and z '■ — 



oo TT 



In the scaling limit when ^ — > oo, we may approximate 



Thus 



TT J-oo 

Thus we get the following scaling form for the real part: 



x'(g,^)~r+'^/(ge,i^) (8) 



with f{x,z) = ^ Jl^dyf{x,y). In particular, the static susceptibility 

x'(g,0)~r+'^/(ge,-oo) 



Thus we identify k = 2 — 7] — if). Note that similar scahng forms apply for the dynamics near 
the classical transition as well, although with a different value for ip and a different scaling 
function f{x,y). Nevertheless, as we argued earlier the function /(x, — oo) should be the 
same for the quantum and classical transitions. 



IV. EXPANSION IN e = 6 - 

As shown below, the upper critical dimension of this model is 6. It is natural to try an 
expansion in powers of e = 6 — d. This was done long back by Aharony, Gefen and Shapir 
T9| and by Boyanovsky and Cardy ||20|. Here, we will review their results and discuss them 



in the context of the general scaling hypotheses of the previous section. 

The e expansion is done using the continuum action Eqn. ||. It is instructive to set up 
the renormalization group so that t is allowed to flow while keeping the strength of the 



randomness fixed. In particular, the results of Ref. p!9|j20|] show that t flows to zero at the 
critical fixed point. When t = 0, the partition function is determined by the particular 
configuration r) that minimizes the action. Clearly the minimum action configuration 
is static. Thus solving the fixed point theory simply corresponds to finding the static con- 
figuraion (j)stat{x) that minimizes the potential energy terms in the action. Thus the fixed 
point theory is entirely classical. 

First consider the Gaussian theory with u = 0. The correlation functions C{k,uj) = 



{(j){k, uj)){(j){—k, —uj)) and G{k, uo) = {(f){k, uj)(j){—k, —oo)) — {(j){k, uj)){(j){—k, —uj)) can be eas- 
ily calculated. The result is: 

G{k,uj) ^ 



k"^ + iu^ + r 

The critical point is at r = 0. The correlation functions introduced in the previous section 
are seen to scale as 

Ci-) - (9) 



GcM-^.Gi\) (10) 

with the correlation length ^ = Thus, in the Gaussian theory, we have i/g = |, r^g = 0, 
TjQ = —2. Note that for the classical problem also, the Gaussian exponents are u = ^, t] = 0, 
fj = —2. This is a trivial illustration of the general point made in the previous section. The 
dynamic scaling is however conventional in the Gaussian theory. 

Now consider a momentum shell renormalization group transformation on the Gaussian 
action. It is convenient to let t flow and keep the strength A of the random field fixed. The 
flow equations for t and r are readily seen to be 

dt 
dr 

-dl=^' 

Thus even in the Gaussian theory, t flows to at the flxed point which is hence fluctuationless. 
The tree-level flow equation for u at the Gaussian flxed point is just obtained by power- 
counting and is 

Thus interaction effects are irrelevant above 6 spatial dimensions and the Gaussian theory 
gives the true critical behaviour. For d below 6, it is possible to construct an expansion in 



powers of e = 6 — c? [19|. To leading order, the one loop RG equations are: 
dt 

dl 4 

where A is the high- momentum cutoff and Kd = j^^- {Sd is the surface area of a unit 
sphere in d dimensions). Again t flows to at the flxed point. Setting t = in the remaining 
equations it is clear that there is a non-trivial flxed point at r* = —^,u* = 35^7^- The 
flows can be linearized around this flxed point and give, for instance, z/ = | + ^ to flrst 
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order in e. Note that this is the same as for the pure problem in 3 — e dimensions. This is a 
general feature of the e expansion - all the exponents characterizing static critical properties 
in d dimensions are the same as the pure problem in d — 3 dimensions. This result however 
is an artifact of the e expansion and is not true in all dimensions. 

The fact that t flows to to this order means that the new non-Gaussian fixed point is 
also fiuctuationless. In fact the flow equation for t has been shown to be exact to all orders 



in e [0,^. Thus at least within the e expansion, the fixed point is fiuctuationless in any 
dimension. 

It was argued by Boyanovsky and Cardy that all the exponents and scaling functions 
associated with static critical properties of the quantum transition were the same as for 
the classical transition. As we have seen in the previous section, this result is true quite 
generally. However dynamic properties (for instance time-dependant correlation functions) 
were shown to scale differently. Within the e expansion, they found that the dynamic scaling 
is conventional at both transitions with zq = l + (0.0185)e^ + (0.0182)e^ + o(e^) and Zd = "^zq. 
Thus, as in the classical case, the e expansion is qualitatively incorrect to describe many 
aspects of the critical behaviour of the random field quantum system well below the upper 
critical dimension. Nevertheless the e expansion provides useful evidence for the claim that 
the fixed point has no quantum fluctuations. 

V. DISCUSSION 

What we have done in this paper is primarily to update the theory of the quantum 
random field models since the pioneering papers in Ref. [0,^], to take into account sub- 
sequent develoments in the understanding of the classical problem. Our main assumption, 
motivated by the results of Ref. p!9|j20[| , was that the quantum transition in any dimension is 



controlled by a fiuctuationless fixed point that also controls the classical finite temperature 
transition. Some of the consequences of this assumption were then examined. The static 
critical properties of the quantum transition are identical to those of the classical transition. 
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The dynamics (which depends on the irrelevant quantum/thermal fluctuations) is similar, 
though not identical. In particular, we suggested that the dynamic scaling is of the activated 
form, with the length scales depending logarithmically on time scales. We have however not 
provided calculational evidence for this suggestion. Evidence for the activated dynamic scal- 
ing cannot be obtained from the e expansion (or from other perturbative approaches like a 
2 + e expansion). Even for the classical transition, the only available evidence comes from 
numerical calculations ||2^ and agreement with the phenomenology seen in experiments [[18 



Here, we have argued that activated dynamics in the quantum case is natural if it occurs 
in the classical problem. We may therefore regard support for activated dynamics in the 
classical transition as some sort of support for it happening in the quantum transition as 
well. 

There are some other important questions that are left open. For most of the paper, we 
have focused on the critical point. The paramagnetic phase should be interesting to study 
by itself. In particular, is there a finite gap all the way till the critical point, or is there a 
Griffiths-phase with gapless excitations in the vicinity of the critical point? It is possible to 



show for the quantum version of a special model introduced by Grinstein and Mukamel 
2^ in c? = 1 that the spin-autocorelation in imaginary time has a stretched exponential form 
g-cTi_ This implies gapless excitations with an essential singularity at zero energy in the 
local density of states. For more realistic models however, even in ci = 1, the autocorrelation 
presumably decays exponentially. The general situation is unclear. Similar questions can 
also be asked about Griffiths effects for the uniform susceptibility in the paramagnetic phase. 
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FIGURES 




FIG. 1. Schematic phase diagram and renormalization group flows for the transverse field Ising 
model in a random longitudinal field for d > 2 in the temper ature(T), transverse field(/i), strength 
of randomness(A) space. There is a transition from a ferromagnet to a paramagnet as any of these 
three parameters is increased. 
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